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Abstract 

We present a non-perturbative method for the study of the 0(N + Inversion of the linear sigma- 
model. Using boson-mapping techniques, in close analogy to those well-known for fermionic systems, 
we obtain a systematic 1/A-expansion for the Hamiltonian which is symmetry-conserving order by 
order. The leading order for the Hamiltonian is evaluated explicitly and we apply the method to 
7T7r-scattering, in deriving the T-matrix to leading order. 



1 Introduction 



In a recent publication |]]J we have presented a non-perturbative treatment of the pion as a Goldstone 



particle within the linear sigma-model using the Hartree-Fock-Bogoliubov (HFB)-RPA theory. Our ini- 
tial interest in a non-perturbative approach was motivated by the desire to find a pion-pion scattering 
equation which obeys unitarity while being consistent with constraints from chiral symmetry, such as the 
vanishing of the a[] scattering length in the chiral limit, —* 0. This is nontrivial since, even in cases 
where the tree- level scattering amplitude fulfills the chiral symmetry requirements, the iteration through 
various reduction schemes of the Bethe-Salpeter equation, such as time-ordered perturbation theory or 
the Blankenbecler-Sugar reduction breaks the consistency Q. Unitary scattering equations are a neces- 
sity, however, when studying medium modifications of elementary cross sections. Indeed for example a 
7r + -7r~ pair in a hot pion gas may form a resonance close to the 2m n threshold or even a bound state 
as was discussed in S [5L JI . Of course such features cannot be treated in chiral perturbation theory. In 
the present work we will therefore generalize the theory to the two-pion problem with application to tt-tt 
scattering and will indeed find a scattering equation which fulfills the requirements of chiral symmetry. 
In constructing the physical pion we have used in our previous work a RPA operator which mixes 
a quasi-pion with a quasi-pion quasi-sigma pair, obtained from a HFB mean-field calculation, in such a 
way that this operator contains the axial charges Q\ as a special limit. The resulting RPA-eigenmode, at 
zero three-momentum, is then identified with the physical pion and we were able to show that it becomes 
a Goldstone mode in the chiral limit. To fully restore chiral symmetry, which is broken at the HFB 
level, it is important to write down the RPA operator in the selfconsistent HFB basis whose vacuum is 



The z's and w are variational parameters and a + , b + are the original pion and sigma creation operators 



ities which are worthwhile mentioning. The main point is that the selfconsistent HFB basis generates 
a mass for the quasipion even in the chiral limit. Despite the fact that they emerge from a respectable 
Raleigh-Ritz variational principle and that Goldstone's theorem is, as it should be, perfectly restored 
once the RPA fluctuations are included this remains a somewhat disturbing feature. In particular for 
the two-pion problem a serious problem emerges since the asymptotic states cannot be defined easily 
in a way which is consistent which chiral symmetry. The HFB quasiparticles can clearly not be used. 
Even though the two-pion RPA, when solved in a similar fashion as the single-pion RPA, will exhibit a 
spurious mode in the chiral limit, this does not imply that the scattering length will vanish. Interactions 
will persist so as to cancel the quasiparticle mass term, in close analogy to the single-pion case. This is 
obviously unphysical. In addition, with explicit symmetry breaking, the threshold will be at twice the 
quasipion mass and not at twice the physical mass. Because of the obvious advantages of a Raleigh-Ritz 
variational theory it maybe worthwhile for the future to consider modifications of the trial wave function 



characterized by the coherent (squeezed) state 
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to guarantee that Goldstone bosons remain as such. In this paper we will not pursue this option despite 
the fact that the approach of |Q] might probably be generalized to the two-pion problem. Instead we will 
identify the origin of the mass generation in the HFB mean field by studying the 0(N + 1) version of 
the linear sigma-model. This will lead to a systematic 1/A-expansion technique which has the advantage 
of being non-perturbative as well (it will also lead to RPA type equations). It will be shown that, to 
leading order, the pion is a Goldstone boson opening the possibility to also derive a ir — ir scattering 
equation to leading order which is physically meaningful. The disadvantage is that the Raleigh-Ritz 
variational principle is lost. In regards to earlier work on non-perturbative approaches to the $ 4 field 
theories in general g § and the l/N -expansion in particular [|To| , p~T| , [l2| we here will however present a 



novel approach. It employs Boson expansion methods and relies on work in nuclear physics where such 
techniques, in connection with the restoration of spontaneously broken symmetries, have been studied 



extensively [13, 14 1. Boson expansion methods for boson pairs in interacting Bose systems have first been 
introduced by Curutchet, Dukelsky, Dussel, and Fendrik [^] and further elaborated by Bijker, Pittel and 
Dukelsky in fl6[| , in close analogy to the techniques known from the fermion case [||, f?J. 
Since the subject of restoration of broken symmetries is full of subtleties, it seems worthwhile to present 
this approach in some detail. Our aim will be to derive a covariant scattering equation for two pions 
to-leading order without taking the large- A limit. This will be achieved by a systematic 1/A-expansion 
of the Hamiltonian of the linear sigma-model. Of course in reality N = 3 and it will be necessary to 
work out corrections to the lowest order tt — tt scattering equation for quantitative studies. Within the 
bosonization framework this can be done systematically and we will leave the evaluation of the l/N 
corrections for future work. 

In detail our paper is organized as follows. First we recapitulate the most important results of our previ- 
ous work, here applied, however, to the 0(N + Inversion of the linear sigma-model. As mentioned above, 
this mainly serves to elucidate the origin of the mass generation at the HFB level. In sect. 3 the boson 
mapping techniques are introduced which lead to the l/N expansion of the sigma-model Hamiltonian. 
The leading-order results are then represented. For the two-pion sector the T-matrix will be extracted 
in sections 4 and 5 where we also comment on the symmetry properties. Conclusions and an outlook 
will be given in sect. 6. 

2 HFB-RPA for the 0(N + 1) Sigma-Model 

As mentioned in the introduction this section serves to identify the origin of the mass generation in the 
HFB-RPA approach of ref. [|J from the point of view of a 1/A-expansion. We therefore consider the 
0(N + 1) version of the linear sigma-model. The difference to the SU{2) x SU (2) ~ 0(4) case lies in the 
fact that one now has an A-component isovector pion field. The Lagrangian density then reads 



C 



2 r 



A 



2 r 



(d.rrf + (d,af\ - ^ [tt 2 + <r 2 j - ^ [vr 2 + a 2 j " + VNca (1) 
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where Ao represents the bare coupling constant, (j,q the bare mass parameter and ir and a denote the 
bare pion and sigma fields, respectively. Chiral symmetry is explicitly broken (in the PCAC sense) by 
the last term in the Lagrangian, ca. 

It is now convenient to define the field operators in terms of creation and annihilation operators as 

7? (x) = / (a q e 1 ^ + a^e'^) 

*( x ) = / i = (b q e^ + b+e-^) (2) 
where the frequency u q , common to both fields, is given by 

= y/l4 + Q 2 - (3) 

Next a canonical transformation is performed for the pion- as well as the sigma- fields by introducing a 
new set of creation and annihilation operators through a Bogoliubov rotation 



0q = X q b q ~ Vqb-q ~ W* 5, 



qO 



(4) 



with u q , Vg, x q and y q being even functions of their argument, and w a c-number. The additional 'shift' 
in the second equation accounts for the macroscopic condensate 



(K) + (bo) _ (x + yo)(w + w*) 



v / (2vr) 3 2 /[ / 



A /(27r) a 2/i 



Ns 



(5) 



To render the transformations canonical the Bogoliubov factors have to obey the constraints 



2 2 i 
X q ~ Vq = 1- 



(6) 



The quasiparticle vacuum |$) ( <3|$) = = ) is now given by the following coherent state 



|$) = exp 



£>i(0)e£ • ^ g + 2; 2 (g)6+6i g ) + ^6+ 



|0>. 



(7) 



were |0) denotes the vacuum for the original basis (a q \0) = b q \0) = ) and z\ = zi = 

It is straightforward to write the Hamiltonian in the quasiparticle basis. The amplitudes u, v, x, y as well 
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as the condensate s are determined by minimizing the vacuum expectation value ($>\H\$>) / {$> | <E»). This 
is equivalent to demanding that the single-particle part of H is diagonal and leads to the well-known 
BCS gap equations and an equation which determines s. These can be cast in a form which contains 
physical quantities: 



£': 



t4 + \ 2 



N + 2 
N 



In + Jjlo + S Z 



3 2 

In + J* + * 



(8) 



where £ n and £ a denote the masses of the quasi-pion and quasi-sigma respectively and I n and l a are loop 
integrals 

wiX , 1 . , i„= f£«, , . (9 ) 



(2vr) 4 qt-El + i^ 



(2tt) 4 q 2 -£l + i^ 



For Af = 3 these results coincide with those of ref. ]!]]. By rewriting the equations for the quasiparticle 

masses as 



£ + ?M r r _ J 1 

^ + 2A 2 S 2 . 



(10) 



one sees that, in the chiral limit (c — ► 0), the quasipion mass remains massive due to the nonvanishing 
difference I n — I a . At first glance this seems to be a 1/N effect, however due to the selfconsistency one 
actually has a power series in 1/N\ 

To render the pion massless one needs to go further and include a set of well-defined RPA fluctuations. 
The procedure is described in detail in ref. [|IJ and thus we will be very brief here. One introduces an 
excitation operator with the quantum numbers of the pion which is, at most, bilinear in the creation and 
destruction operators 



Qt = Xla+ - Y^a + £ [x 2 n (q)^at q - Y*(q)^ 



qQq 



and similarly for the sigma field 



Y*(q)(3- q (3 q ] + £ 



—q 



at. 



(12) 
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The RPA ground state is determined by Q n \RPA) = Q a \RPA) = 0. Within the equation of motion 
method ]?], [lTj the RPA frequencies can now be determined which finally leads to the following equation 
for the physical pion mass 



2_£ , 2Ag [g - g] [£^(0) - S^K)] 

7T ~ + »t ox, 2 

^ 1 - ^S w (m2) 



(13) 



where ^^(p 2 ) is the quasipion-quasisigma bubble given by 



£trt(p 2 



d A q 



(2tt) 4 - £2 + irj (p _ 9 )2 _ 52 +iT/ - 



(14) 



In the chiral limit (c = 0) the Goldstone theorem is now manifest since eq. fll3|) has a zero-energy solution. 
Here again one sees that the RPA has brought about a power series in 1/N which, in the chiral limit, 
cancels exactly the mean field contributions in such a way that the symmetry is restored. 
After the inclusion of RPA fluctuations one obtains for the sigma mass 



2\ A s 2 



iA^S 7r7r (m 2 7 )S (T(T (m 2 r ) 



(15) 



with the 7T7T and aa bubbles given by 



iY, aa {p 2 ) 



d A q 



(27r) 4 q 2 — E 2 + irj (p — q) 2 — E 2 + irj 
d A q 1 1 



(2tt) 4 q 2 -£ 2 + irj (p - q) 2 - £ 2 + i V 



(16) 



It should be mentioned that the condensate s does not receive any contributions from RPA fluctuations 
in this approximation. 

In spite of the fact that the Goldstone theorem is fulfilled by the formalism described above it has 
a severe draw back for the it — it scattering problem. As mentioned in the introduction there will be 
interactions in the chiral limit due to the fact that the quasipions remain massive. 

It is now interesting to consider the limit (N — > 00). By inspection of eqs. (^),(|l^)and ( |l5|) the HFB-RPA 
approach yields the following results 



£ 2 = - 
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- s +2\ls 2 



m 2 - f 2 + 2A 4 -7 2 ^{ml) 

where the Ti nn bubble is built out of the asymptotic states 

ZUP 2 ) = ~ij ^Ar(?)Ar(p - q) (18) 

which are 'acceptable' from the point of view of chiral symmetry. Obviously the pion is now a Goldstone 
boson already at the mean-field level. The RPA fluctuations do not survive in the large-A^ limit in the 
case of the pion while for the sigma they are still partly present. In fact, as will become clear in sect. 3, 
the result obtained for the pion mass corresponds exactly to the Hartree-Bogoliubov solution since terms 
induced by the Bose statistics disappear. The above results have been obtained previously by several 



authors using various methods |T^, 11]. What we have presented so far is yet another way of deriving 
these results. They merely serve as a motivation for introducing a systematic 1/N expansion in the next 
section. 

First we wish to address another issue, however. In contrast to the HBF-RPA approach in 0(4) the 
large- N limit permits a relatively straightforward regularization. The bare coupling constant Ao and bare 
mass are replaced by renormalized quantities: 

a2 = i $ T m f = * + a2jLo(a)) + a2 ^ o(a) (19) 

where Lq and Kq represent a logarithmic divergence and a quadratic divergence, respectively, and A is a 
four-dimensional cut-off. Defining the regularized pionic tadpole and the irir bubble as 

h = I n - K (A) - mlLo(A) t^{p 2 ) = ^ n (p 2 ) - L (A) (20) 

the physical pion- and sigma masses are rendered finite 



2 c 2 , \2 

m% = - = fi z + A z 



U + S 2 



2AV 

s ' [1 - A 2 S^(m2)] • 



ml = - + t- 7"- , (21) 



In addition, the well-known Ward identity for the anir vertex [18| is obeyed: 



> > 2 A s 

m a ~ m n = T - , s. 22 

[1 - A i! S 7r7r (m z )J 
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A comment is in order here. From eq. ( |19| ) one can easily see that the renormalized coupling constant 
A vanishes when the cut-off A is removed and hence the theory becomes trivial. The triviality of the 
A<£ 4 -theory in the 1/N approach has been thoroughly studied (see for instance [11, ^]) and, to leading 



order, it has been shown in [|ll]] that the vacuum has the same group of invariance as the Lagrangian. The 
same is evidently true here. This problem may be avoided, however, by assigning a physical significance 
to a finite cut-off A, arguing that the sigma-model only describes physics up to a certain scale. In this 



spirit it was shown in ref. [18] that, with a cutoff version of the A$ 4 -theory, the nontrivial vacuum can be 
stabilized. Another interesting possibility is the introduction of more degrees of freedom (vector bosons 
in the case of A<3? 4 [^(J) which may have its own difficulties, however. Clearly more work will be needed 
on this point. 

3 Dynamical Boson Mapping and the 1/iV-Expansion 

In this section we present an approach to the 1/A^-expansion in the cut-off version of the 0(N + 1) 
sigma-model inspired by well-known boson mapping techniques used in nuclear physics || ^, [l4[]. This 
method is systematic in the sense that it represents a well-defined expansion in powers of l/\fN and 
no a posteriori limit N —* oo has to be taken. In fact, the RPA is based on the so-called quasi-boson 
approximation which, in the nuclear physics context, means that Fermion pairs are replaced by ideal 
bosons. Analogously the RPA for bosons, as applied in the last section, implies the bosonization of boson 
pairs. Boson expansion techniques for interacting bose systems have first been studied by Curutchet, 
Dukelsky, Dussel, and Fendrik (CDDF) in |1§ and by Bijker, Pittel and Dukelsky (BPD) in @. We will 
apply this technique to obtain a systematic 1/N expansion of the Hamiltonian of the linear sigma-model. 
It has been discussed at length by Marshalek et al. [|b|, 14 1 how such an expansion preserves at the same 



time the symmetries order by order. We are specifically interested in the it — n scattering problem for 
which we want to establish a Lippmann-Schwinger type of equation which fulfills unitarity as well as 
constraints from chiral symmetry such as the vanishing of the it — n scattering length once the explicit 
symmetry breaking is removed (m n — > 0). To proceed one introduces ideal bose operators (A+ ,A qtP ) 
and bosonizes pion pair operators as first laid out by Holstein and Primakoff |^l| (in nuclear physics it 
has become known also as the Belyaev-Zelevinsky method [|22}[). The method consists of the following 
mapping 



a q a p = {a^a+^j 

a+a p = (A + A) qp (23) 
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where the new operators (A+ p , A q ^ p ) are ideal bosons and obey the usual Heisenberg-Weyl algebra: 



A A' 



p,q 



A + A + 

■™-m,n > p,q 



{.8 m ,q8n,p 8 n ,q8m,p) 
\A m ,n i A p ^ q \ — 0. 



(24) 



With respect to the more familiar Fermion case one should notice the change of sign on the rhs of the first 
equation in eq. (|24D and under the square root in eq. (p3|). In addition, the operators A q ^ p are symmetric 
under exchange of indices 



A —A 



A + = A + 

q,p Vil 



(25) 



rather than antisymmetric as in the case of Fermion bosonization. The approach of Belyaev-Zelevinski 
to the bosonic mapping requires the realization of the original algebra for pairs of pions 



0*10*2 , 


a£a£ 


a \d 2 


a 4 




, a^at 



N ((5l3#24 + 5l4<523) + #1304 0,2 + £230*4 0*1 + 0140*3" 0*2 + £240*3 0*1 

0130*20*4 + £230104 
[0*10*2 , o 3 a 4 ] = 0. 



(26) 



The reader may verify that this is in fact true order by order, using the Holstein-Primakoff mapping in 
eq. (H). 

To write the Hamiltonian in the (CDDF-BPD)-representation we need to define the sigma field. Since 
the latter is to be constructed perturbatively, it is legitimate to use a form similar to eq. (pi) 



a x 



d 3 i 



(&q 



e iqx + 6+e" iqx 



(27) 



with the sigma-quasiparticle mass E a to be defined later. With these definitions the Hamiltonian takes 
the following form: 



H 



n + 



(Bl{q) + u> 2 q + \lh) btb q - l -(e 2 M) 



u? q - A§/o) (bpt q + b- q b q 



+ Jd 3 q e (q){A+A) qq + 



Ao(g) 



+ 



Jdx ^(- 2 (x) :+d*(x)) 2 -VN, 



(A + VN + A+a) + (y/N + A+AA 
Nca(x 



-q,q 



(28) 



8 



where the colon indicates normal ordering and the various functions are given by 



Ho 

eo(<?) 



NX 2 j2 



with Iq 



A (g), 



d 3 q 1 

(2vr) 3 %T q 

Ao(fl) = ^U>. 



(29) 



The normal-ordered pion fields can be explicitly written as 
d 3 giA 2 e i(qi+q2)x 



: 7r x : : 



(27r)V4^i^2 



(30) 



Clearly the Hamiltonian is an infinite series in powers of the new bosons as a direct result of the Holstein- 
Primakoff mapping. It then takes the form 



H = + WW + «(2) + WW + + 



(31) 



The superscript in each term designates the power in the boson operators and these are either the sigma 
bosons b q ,b q or the newly introduced pion-pair bosons A, A + . In fact, in addition to this expansion 
in the power of bosons there exists an underlying expansion in the square root of the number of pion 
charges, yN. The latter is particularly suitable for a perturbative treatment of the dynamics generated 



by the Hamiltonian eq. (28). The fact that N is arbitrary insures the preservation of the symmetries of 
H at each order. The perturbative calculation of all higher orders is performed in a selfconsistent basis 
which we will establish first. Since there are two bosons which are able to condense in the vacuum, the 
natural choice for the trial wave function of the ground state is a coherent state of the form 



exp 



Y,d q A+_ q + (b )b+ 



|0>. 



(32) 



To fix the parameters d q and (bo) one needs to apply the Ritz variational principle to the expectation 
value Tt° = ^^y-- The latter can be evaluated and gives 



(2it) 3 N 



d 3 q d 2 



(2tt)3 q N 



-^r- (I + K ) + - (Jo + Ko) +^2- + ^T~^ 



(33) 

The parameter (a) = ^ 2 + )3 2£ ' aS d enne d * n e( l" ©' * s n °thi n g but the expectation value of the sigma 
field in the coherent state while Kq is given by 
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d 3 q { d q + y/N + q ] :v 



N J (2vr) 3 



2cj„ 



(34) 



The minimum of H determines the values of d q and s = (a) /v N, as solutions of the following equations 



e n (q) [2d q ^N + <P q ) + A n (q) [2d q + N) = 

+ 2uj q A n (q) = - 



(35) 



where A n (q) and e 7r (g) are given by 



2w„ 



I + K + s z 



(36) 



The first equation in (35) can be reexpressed in the more suggestive form 



AT 



(37) 



From the canonical form of this 'gap equation' one can deduce the pion quasiparticle energy £ n (q) as 



£ n (q) = Jel(q) ~ A>(fl) = Jtf + + ^ [I + K + ^ 



(38) 



One can also deduce from eq. (36) that the shift d q for the pion pairs can have a simple scaling, namely 
as v/iV, which shows the physical similarity of this parameter to the scalar field condensate (a). In what 
follows we will use a rescaled shift parameter 



d q = dg/VN 



An interesting identity which allows to transmit the selfconsistency of the gap equation in eq. ( pq ) to the 
pion quasiparticle mass reads 



M?) = ( d i - \A + d l 



(39) 



from which one can finally extract the Hartree-Bogoliubov (HB) mass of the pion, 



m n = £ n (0), 
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and the value of the condensate via the two coupled equations 



m n = ^0 + A 
C 

s 



(40) 



where !„■ = Jo + Kq denotes the quasi-pion tadpole. 

Having fixed the quasiparticle basis we can now expand the Hamiltonian in eq. (p 
shifted bosons operators A qiP and q defined as 



in powers of the 



A -A 

9iP — ItP 



with 



Nd q 5(q+p) f3 q = b q - (b ) 5(g) 

= W) = o. 



(41) 



Knowing the scaling of all parameters in the Hamiltonian the latter can therefore be expanded without 
ambiguity according to 



H = NH<® + VNH^ + + ^=H^ + -H^ + .. 

Vn n 



(42) 



where each order is obtained by expanding the square roots in eq. (|23|). Using the parameter differenti- 
ation techniques of operators (see appendix 1) one can write down the contributions to the Hamiltonian 
to leading orders in 1/N, namely and H^: 



'(27T) 3 

~2^r 



^qIttS + AgS^ + /igS 



+ 



d 3 q 



LUgdq + 



(fa + Pi 

A„{q) (d q + yjl + df) 
2jl + d\ 



(A+_ q + (43) 



# (2) = J d 3 q£ a (q)P+/3 q + J d 3 q£ n (q) (l + dfj (I+i) 



+ I d 3 qd 3 p 



M«) V 1 + d l ~ V 1 + d l f( dA + A+d) q p A m + At p (dA + AU) p ( 



d\ dp 



+ 



d q A n (q) / 1 



l + d? q - Jl + dl 



d\ ~ d 2 p Wll + dl 



d? 



\ (dA + A + d) (dA + A+d) 

di \ 'q,p\ 'p,> 
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d 3 qid 3 q 2 d 3 q 3 



+ f 



^(27T) 3 8£ a (q 3 )io( qi )u;(q 2 ) 
f d 3 qid 3 q 2 d 3 q 3 d 3 q 4 5 3 (qi + q 2 + q 3 + Q4 



1.2 



(27r) 3 v /16w(gi)w(g 2 )w(Q , 3)^(Q , 4) 



<r(gi + Q2 + g 3 )ri, 2 /% + 0I 3 \ [A 

Al, 2 + iii_ 2 



1,-2 



Tl,2 T3 4 



^3.4 + a: 



3,-4 



(44) 



with the following definitions for the sigma quasiparticle mass E a and for the factors r$ ,•: 



ml + 2A^s 2 



di + a/1 + e£ 



2(di - d*) 



(45) 



The gap equation which we have derived in eq. (^) eliminates the linear part of the Hamiltonian, i.e 
ii^ 1 ) = thus defining the HB basis. The leading contribution in the expansion is therefore of order 
one. As we have seen, an appropriate mapping of the boson pair operators into ideal bose operators 
automatically produces the correct result. 

The Hamiltonian above can be recast in a more compact form which is particularly suitable for the 
forthcoming considerations. To achieve this a set of new operators B qp , B q)P is defined which corresponds 
to a rotation of the set of A+ p , A q>p operators, according to 



(46) 



In order for the B+ p ,B qtP operators to obey the same algebra as the A+ p ,A qjP , the coefficients G qtP and 



H q ^p have to be chosen such that 



G 



q.p 



H. 



q.p 



1. 



(47) 



Explicitly the functions G qp and H q)P are given by 



dqyjl + dj ± dpyjl + d\ 

d q ± d p 



(48) 



With these definitions one verifies that the symmetry of the operators A£ , A q>p with respect to the 
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interchange of the indices is also obeyed by the new bosons B£ , B qp , since G q ^ p and H qp are symmetric 
under the interchange of the indices. The quadratic part of H now reads 



H 



(2) 



d 3 q £ a {q) /3+[3 q + / d 3 qd 3 p £ v (q) B qp B q>p 



d 3 q 1 d 3 q 2 d 3 q 3 S 3 (q 1 +q 2 + q 3 ) 



fa + /£ 3 B lt2 + B+ 



+ f 



f d 3 q 1 d 3 q 2 d 3 q 3 d 3 q 4 5 3 (qi + q 2 + q 3 + q 4 ) 



(2ir) 3 v / l6£ n (qi)£ n (q 2 )£ n (q 3 )£ 7T (q4:) 



B 1>2 + B+ x 



B3.4 + B_ 



3,-4 



(49) 



The next task is to diagonalize each order of the Hamiltonian. Parts with odd powers in the operators 
have to be eliminated by this procedure, as was the case for iJW . Expressing in diagonal form leads 
to an expression with uncoupled oscillators as will be seen below. 

The diagonalization of can be performed by using a general Bogoliubov rotation which mixes the 
configuration of the quasisigma Boson and the bosonized pair of pions according to 

Qt (p) = uP (p)/3i - (p)P-p + ^= E K 2) (8 9) B tf-? ~ v ? ] ® fl S -*H**] • (50) 

The operator can be considered as an RPA excitation operator. The rotation in the RPA operator 
in eq. ( [50] ) is performed under the kinematical constraint 5(qi + q 2 — p), where p denotes the total mo- 
mentum. 

Using Rowe's equation of motion |[7]] 

(RPA\ \5Q u (p) , \H^ , Q+(p)ll \RPA) = n u (p)(RPA\ [5Q u (p) , Qt{p)} \RPA) (51) 



and keeping in eq. (51) the full RPA ground state defined by Q V \RPA) = would lead to a selfconsistent 
form of the RPA equations (see |23| ) if higher than second-order terms are included in the expansion of 
H. Retaining only linearizes the equation of motion and one obtains the familiar form of the RPA 
equations 

d% ) ( W>&) \ = nMN ( KM) , (52) 

Bpiqi,q2) Apiqucfit) J \ V u {p,q 2 ) J y V v (p,qi) 

where Ap, Bp are 2x2 matrices given by 

Ap(qi,q 2 )=\ ^ / r „ ^ _ _ N1 1^(91-92) + %(<fi,<f 2 ) (53) 

[£-K{qi) + £tt{P - qi) 



13 



with 

y/2Ags A| 

V \/ (2Tr) 3 S£^{p)£ 7v (qi)S 7v {p-qi) (2tt) 3 ^/ 16£,r (<?i )£,r (<? 2 )£■* (p-Qi )£■* (p-<?2 ) 



%(gi,g2) 



(54) 



respectively. The 4x4 norm matrix as well as the two-component vectors Z^, V v are given by 



V = | " ° 7,= 1 ° 

-h j \ 1 

One can now solve the eigenvalue problem and extract the RPA frequencies to obtain 

Km - e * + i - Agiwns®) P ( } 

with Ti nir {Vt^{p)) being the Lorentz invariant tttt self energy of the sigma given by 

V (0 2 (^\ f d ' iq £ ^ + S ^P~l) I (K 7 \ 

WKW) ~ J (27r)3 2S n { q )E K (p-q) Q 2 (p) - (£ n (q) + ^(p-,)) 2 ' ^ 

It should be noted that, in contrast to (|l^), where the quasipion was massive in the chiral limit, here 
the E,,-,,- self energy in (|57]) is built on massless Goldstone pions. Using the usual orthonormalisation 
condition of the RPA states, i.e {RPA\Q U Q+,\RPA) = S uu i, one obtains the following solution for the 
RPA eigenvectors 



UP(p) 



(2tt)2 K 



^2 E(p,n u )^/2£J^jn u -£ a (p) 
(2tt)I 1 1 

1 ^7T7r-^7r7r 1 1 



2 H(p, «„) ^A£ n {q)£ n {p-q) Q„ - £ n {q) - £ w (p- q) 

1 VtTTT — J-7T7T 1 1 



2 H(p, a,) y/4£ n {q)£ n (p- q) fi„ + ^(0) + ^(p - g) 



(58) 



where V^^cr and V^tt—^^^tdp) stand for the 7T7T — * a and 7T7t — * irir tree level transition matrix respec- 
tively. These as well as the function S are given by 
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(2vf 2 
Ml- 



i 



2 5S vr7r (^,p) 



1/2 



90, 



(59) 



This concludes the discussion of the leading order dynamics. The physical observables, at this order, are 



given by the expressions in eq. (40) for the pion mass and the condensate and by eq. (|56| ) for the sigma 
mass. These were the same results obtained at the end of the previous section. 

In the next section some remaining subtleties in relation with the properties of the RPA approximation 
will be pointed out. 

4 The 7r — 7r Goldstone Mode and its Dynamics, Sum Rules. 

In the chiral limit the RPA solutions have to yield a spurious mode corresponding physically to the case 
of two non-interacting Goldstone pions. This mode has to occur if the massless pions have zero total and 
zero relative momentum. Such a solution is not manifest from expression (|56|) which has been obtained 
by solving the RPA equations in the sigma channel via Feshbach projection. In the chiral limit and for 
vanishing total momentum p the RPA frequency will just become the sigma mass as to be expected and 



as can also be seen by comparing eq. (56) with eq. (20). That a Goldstone mode (i.e. a Goldstone mode 
for the two pions) must also exist, can be deduced from the following symmetry considerations (see also 
section 5). 

Since, without explicit symmetry breaking, the sigma-model is invariant under chiral transformations, 
the Hamiltonian must fulfill the following commutation relations 

H,Ql] = [H,Q 5 Q 5 ]=0 (60) 

since Q\, the i-th axial charge, is a generator of the symmetry. The second relation is most interesting in 
the present context since it involves two-pion degrees of freedom. As for single Q\ operators, the product 
of two axial charges can also be mapped into the bosons A + ,A or equivalently B + ,B. This allows to 
express S = Q5Q5 in an infinite series in powers of the new bosons 

S = N 2 S ( V + N 3 / 2 S^ + NS { V + >/NS® + + -L + ... (61) 

y/N 



The S(p) are organized according to the power p of the bosons (3, j3 + , B, B + . It is important to restate 
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that, as in the case of the Hamiltonian expansion in eqs. (42.43,44), no normal ordering is assumed. 
Therefore a given term S^) may very well contain p — l,p — 2, .. powers of the bosons. The real expansion 
parameter in eq. ( f42|) as well as in eq. ( |6"l"l) is y/~N. 
The commutation relation eq. (|60|) will thus read 



[H, S] =^^3-1^-1) 



(62) 



According to the basic observation that the polynomial is zero if and only if each of the monoms is van- 
ishing, one can thus extract less stringent commutation relations between parts of the total Hamiltonian 
H and the symmetry operator S. These are simply 



C (P) = \ H (1) 5 (p+l)l + \ H {2) 5 (p)l + ....+ \ H (P+1) j 5 (1)1 = 



Vp > 0. 



(63) 



The relevant commutation relation to the order we are working is given by CO) = and reads explicitly 



ffW , 5< 2 >1 + , so)} = 



(64) 



Recalling the fact that HO) is eliminated by the construction of the variational (HB) basis, one is left 
only with 



(65) 



This proves that a zero mode must exist. For the discussion below it is useful to give an explicit form 
of the symmetry operator. Since S^ 1 ) contains only unit powers of the bosons j3, /3 + , B, B + , the RPA 
operator Eq. (50) is therefore the most general ansatz with single boson excitations containing so-). 
Explicitly, SO) can be constructed as 



with 



U sp (0,0) = -V sp (0,0) 



+ 



(27r)a a (2f ff (0))-3f ff (0) 
-(2vr) 3 4^(0) 



(66) 



(67) 



Physically this operator creates a zero-momentum sigma meson and a pion pair of vanishing total- as 
well as relative momentum. 

As mentioned above, it is instructive to express in the RPA basis. We shall follow the usual 
steps that can be found in the literature (see for instance M). First it is easy to see that, by introducing 
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the RPA matrix, can be written in matrix form as 



#(2) = _I Tr f A _ B ) + - f/3 + -^=B + -^=B 



[a b\ 



B A 



t » \ 



7t B 



Following Q we introduce a matrix T and denote the RPA matrix by 5: 



T 



U V* 
V W 



A B 
B A 



(68) 



(69) 



Then the closure relation of the RPA basis \M reads 



TAAT+ = AT 



(70) 



where AT is a two by two norme matrice defined by substituting 1 to Id in N of eq.(|55|). This allows to 
reexpress the original boson operators in terms of RPA excitation operators as 



B + 

q,p 



B, 



q.p 



V2j2(vP*(-q,-p)Qt(-p-q) + UP(q,p)Q v (p + . 



(71) 



u>0 



Now, using the orthonormalisation condition of the RPA basis which in compact notation reads 



T+A7T = M 



(72) 



one can extract, with the help of the RPA equation, the components of the eigenvectors U u (p,q) and 
Vv{p,q) as defined in eq. (|55|) and explicitly given in eq. (|58|). Finally, with the RPA excitation operators 
as defined in eq. (IbO) the Hamiltonian takes the following form 



H m = -\ T r{A-B) + ~ (Q+ Q)T+5T , 

2 2 \ Q- 



(73) 



In the exact chiral limit we have seen that a zero mode (with Q sp = 0) is present in the set of RPA 
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solutions which requires special treatment when inverting the RPA. The procedure is standard and can 
be found in textbooks (see for instance J?]]). Taking into account the zero mode which appears in the 
chiral limit, the completeness of the RPA basis is again recovered and the Hamiltonian can be expressed 
as follows 



i i (s (1 T 

--Tr(A-B) + ^ E ^(p) + E ^(P)Qt(p)Qu(p) + ^f- 

p, v>0 p, u>0 



(74) 



is the 0(1) symmetry generator given in eq. (61) and 2q is the moment of inertia per unit volume 
(with V the volume of the system) which can be calculated using the Valatin-Thouless equations [24]: 



#(2) )C (D] = S (D 



(75) 



The operator C^- 1 ' is an anti-Hermitian operator proportional to the canonical variable conjugate 
to the zero mode induced by such that 



C« = UnT^ 



with 



(76) 



The Valatin-Thouless equations in eq. ( |75| ) are sufficient to determine both and the moment of 
Inertia Zq/ unit volume. An explicit calculation gives 



T (i) 



(27r) 3 2s 2 m* + 1 
c2 



B, 



oo 



V(2^) a 2g CT 

(2^)328 777^ 



(/5o + " 0o) 



-o = -r 



(77) 



The expression for the moment of inertia suggests the very intuitive picture that the inertia of the vacuum 
against a chiral rotation is simply given by the amount of the spontaneous breaking via the condensate, 
s. Of course in general the volume V in eq. ( |74"|) is infinite and thus the kinetic energy term of chiral 
rotation in ( |74| ) is zero. However, in relativistic heavy ion collisions, situations may occur where there 
are small droplets of pionic matter with spontaneously broken chiral symmetry, surrounded by regions 
in the symmetry restored phase. In such cases the inertia is finite and we must keep the corresponding 
term. 

Before closing this section, we wish to address briefly the question of the energy-weighted sum rule 
(EWSR). For any hermitian operator, F, this sum rule, to leading order in the 1/iV-expansion, is given 
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by 



m x F = Y,^MF\RPA)\ 2 = l(nF,[H (2 \Fm (78) 

where l^) is the coherent state specified in eq. (^). The proof goes along the standard arguments which 
can be found for instance in J7| 3 |25| . 

Among other things such sum rules can serve as a test for the correctness of numerical calculations. For 
example, a particular simple result is obtained with F{x) = <r(x) for which the EWSR gives rrip(q) = 1, 
where q is the three-momentum carried by the a-field. Further sum rules with different excitation 
operators will be discussed in a forthcoming paper. 

This concludes the formal discussion of the 1/N expansion and the properties of the ensuing RPA 
problem. We will now turn to the construction of the tttt T-matrix. 



5 The 7T7r-scattering equation 



A 7T7r-scattering equation can be deduced from eqs. (51, |52| ) by eliminating the sigma subspace with the 



help of a Feshbach projection. This procedure adds to the original 7T7r-contact term of the Lagrangian 
an effective part which corresponds to s-channel sigma exchange. To leading order in the l/A r -expansion 
the projected RPA equation is given by 

d 3 q\([S n (qi) + £ n (p- qi)]5(qi - q 2 )I d + M^u,P,qi,q2)) \ \ 2) / x = & v {p) 

\ Vv {p,qi) / 




Here Id is the two by two identity matrix and the interaction matrix A is given by 



k(£l v ,p,qi,q 2 ) 
d{p,qi,q2) 



A 2^ 



£ 2 M 



d(p,qi,q2) 



Slip) \1 1 J (2^) 3 

[I6£n(qi)£n(q2)£n(p- qi)£-n{p- <72)n 1/2 



(80) 



From eqs. (f7|) and ( p0| ) it is now straightforward to construct the corresponding two-pion propagator in 
the scalar-isoscalar channel: 

G^{t-t',x x ,x 2 ,xz,x±) = -i( J RP^|T(vf(f 1 )^(x 2 )) i (7f(x 4 )7f(£ 3 )) t , \RPA) (81) 



where T denotes the time-ordering operator and \RPA) is the RPA-correlated vacuum. The Green's 
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function only depends on the time difference t — t' but, in spite of being non-covariant, still contains the 
full dynamics at this order since t or u-channels only appear in next order. Taking the Fourier trans- 
form of G T7T and keeping the relevant part to leading order of the expansion in the new bosons B, B + gives 

G^(t-t',q 1 ,q 2 ,q3,q4)=Nd{qi,q 2 ,q 3 ,q i )(RPA\T(Bg u g 2 + £± gi ,_ g2 ) t + B± qs ^ ^ \RPA) 

(82) 

The Fourier transform of the time difference t — t' defines the center of mass energy E . At this order 
the two-pion Green's function is related in the usual way to the T-matrix 



G-rir (E,qi,q 2 ,q 3 , q&) = Gl n (E, q x , q 2 , q 3 , q&) + Gl n (E,q\,q 2 )T(E,p)G° n (E,q 3 ,q\) 
and a straightforward reduction then gives the T-matrix as 



(83) 



dkidk 2 dk 3 dk4 G° 7rn 1 (E,q 1 ,q 2 ,k 1 ,k 2 ) <SW(£, 9i, <?2, 23, £0 G" n l (E,k 3 ,k4,q 3 , q\) 



(2n) 



hpS(qi + q 2 -q 3 - q\)T^ n (E,p) 



(84) 



where p is the total three-momentum of the pion pair, p = q\+ q 2 = q 3 + q 4 - On the other hand G"^ 7 
the free two-pion Green's function defined as 



is 



Glnit- t' ,q\,q 2 , q 3 ,q4) = N d(qi,q 2 ,q 3 ,q 4 )(ip\T [B, 



>< n .<i2 + B_ qi _ q2 ) t yBq 3>q4 + B_ q3 _ q4 j^ 



G°^(t - t', q h q 2 ) [S(qi - q 3 )6(q 2 - q 4 ) + 6(qi - q*)5(q 2 - q 3 )\ (85) 



where is the coherent state defined earlier. In energy space one has 



G® w (E,q\,q 2 ) 



2£ K fa)S v fa) E* - (S n (q\) + £^q 2 )f + i V 



(86) 



Using the Heisenberg picture, the reduction in eq. ( p4|) can be performed explicitly with the help of the 
RPA representation in eqs. (71, 74) as well as the RPA solutions in eqs. ([56], 57, ^ , |59| ). This yields the 
following spectral representation for the T-matrix 



T n7T (E,p) = - 



v>0 



V^(VL v ,p) 



i(n„,p) 



i 



E - n v (p) +ir] E + n„(p) - ir) 



(87) 



where the V ww and H functions have been given in eq. (^9|). An equivalent way of solving the scattering 
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problem in the RPA is to extract the bare 7T7r-interaction V nn which serves as the kernel of a Lippmann- 
Schwinger equation: 

T n7T (E,p) = V n7T (E,p) + lj -^V^{E,p)G^{E,p,q)T^{E,p) 



with G® n as given in ( |86|) and the bare interaction V nw obtained by reducing the Green's function in 
eq. ( |8l| ) to lowest order in the dimensionless coupling constant, Xq. The solution of the above equation 
has a simple algebraic form 



T(E,p) 



lV wn (E,p)T,^ 7r (E,p) 



A 2 = N £Z(0)-jt(0) w . th f 2 = Ng 2 m 



where Ti ww (E,p) denotes the 7T7r-bubble defined in eq. (|57|), The two forms of the T-matrix given above 
are equivalent. This can be verified by taking the residues at each pole of T nn from the expression in 



eq. ( |89| ) which then directly leads to the spectral representation in eq. (p7|). From eqs. ( |82| , |8^ , 84, |89| ) 
and from the imaginary part of the total two-pion propagator ImG^n one deduces that the spectrum 
starts at 2m n i.e. it is gapless in the chiral limit, in agreement with the Goldstone Theorem (see sect. 4). 
The a -propagator can also be calculated starting from the related two-point Green's function 

D a (t 1 -t 2 ,x 1 ,x 2 ) = -i {RPA\Ta(xi,ti)a(x 2 ,t2)\RPA) (90) 



which, after Fourier transformation, can be written in spectral representation as 



D a (E,p) = -J2 



u>0 



i 



E - n„(p) +ir) E + n u (p) - ir] 



(91) 



Alternatively one can obtain D a by solving a Dyson equation with the mass operator containing the 
7T7T-RPA fluctuations. This leads to 



D a {E,p) 



E 2 -81 



2Xy ^ n (E,p) 
1 - X 2 ^{E,p) 



-i 



(92) 



and again these two forms of D a are equivalent, as can be seen by taking the residue at each pole in 
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eq. ( |92| ) which yields the spectral form (^). From these observations the scattering matrix T n7r can be 
reexpressed as 

W) = (93) 

where p 2 = E 2 — p 2 . As we have seen above to this order in the 1 /N expansion the pion is obtained in 
the Hartree approximation, and D n is therefore given by 

Ar(p 2 ) = 2 1 2 - (94) 

Thus one pion can be taken soft (i.e. = 0) in which case the above expression reduces to a Ward 
identity, linking the pion four-point function to the sigma and pion two-point functions ||18|| . The expres- 
sion in ((9^) is more general, however, and includes the soft-pion limit as a special case. Finally it can be 
seen that, at the physical threshold E 2 = Am 2 r ,p = 0, the T-matrix is directly proportional to the pion 
mass and hence vanishes in the chiral limit. This is, of course, required by the Goldstone nature of the 
pions. Away from the chiral limit it is interesting to compare the a{] scattering length resulting from the 
T-matrix in eq. ( |89|) with the well-known tree level results. To lowest order in the coupling constant and 
for N = 3, the a[] scattering length can be read off from eq. (|89|) and yields 

9 m n 

a 0tree ~ 32^ \ ' 

while the standard tree-level value is a[] = ^^yf- The difference stems from the fact that, to leading- 
order in 1/N, only s-channel contributions appear while the standard Weinberg tree-level result contains 
t- and u-channel pieces as well. These contributions are known to be repulsive and, in the 1/N expansion, 
they only enter in next-to- leading order. It is tempting to obtain the a{} scattering length for the full 
T-matrix. Unfortunately for any reasonable choice of parameters (the cut-off, A, and the physical sigma 
mass, £ a ) no acceptable fit of the scalar-isoscalar phase shifts can be obtained from eq. (|89|). This is to 
be expected since the RPA resummation of the scattering series will increase a[] further from its tree level 
value, which is already too large. This hints to the fact that higher orders in the l/Af expansion will be 
crucial to get a quantitative description of the empirical s— wave phase shifts. It will be very interesting 
to work these out explicitly, exhibiting u- and t- channel contributions. This will be done in the future. 

6 Summary and Outlook 

In an attempt to obtain a symmetry-conserving 7T7r-scattering equation, consistent with unitarity, we have 
investigated the cut-off version of the 0(N + 1) linear sigma-model. Increasing interest in connection with 
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two-pion correlations in hadronic matter makes such studies particularly timely. On the one hand, as is 
well known, two-pion correlations play a crucial role in the N-N interaction and studies of the in-medium 
corrections are still in their infancy f2,|26]. On the other hand, in purely pionic matter, the 7T7r-interaction, 



due to Bose statistics, may become strongly enhanced such that, in the scalar-isoscalar channel, two-pion 
bound states are likely to appear at higher temperature. Obviously, such features cannot be treated in 
chiral perturbation theory and non-perturbative schemes, respecting constraints from chiral symmetry, 
are called for. The present paper provides a novel approach by applying well-known Boson-expansion 
techniques. That such techniques, which are highly developed in nuclear physics, can be very useful 
for interacting Bose systems was first demonstrated by (CDDF-BPD) [|l^, We have applied this 
approach, for the first time, in the context of field theory. In a systematic 1 /^-expansion it has become 
possible to derive a 7T7r-scattering equation, to lowest order in 1/N. This scattering equation is of RPA 
type and fulfills all required properties such as unitarity, Ward identities and standard RPA sum rules. 
The formalism is manifestly non-covariant. However the full covariant theory collapses to a two times- 
theory in this order of the 1/N expansion. Therefore no breaking of covariance occurs in the final 
results. Furthermore special care had to be taken in dealing with Goldstone modes, which is crucial 
for 7T7r-scattering in the chiral limit. Due to the symmetry-conserving character of the approximations 
these aspects are treated correctly in the present approach. The lowest-order equation disregards u- and 
t-channel exchanges in the I = J = channel and no physically reasonable parameters (cut-off and sigma 
mass) could be found to reproduce the empirical ttit phase shifts. Contributions from u- and t-channel 
exchanges , possibly curing this problem, enter only in next-to-leading order and will be studied in a 
forthcoming publication. 

One can envision several other applications of interest. A direct extension of the present study would 
be the inclusion of fermionic degrees of freedom to treat the tt — N system on a new level, including 
non-perturbative irir rescattering effects systematically. Another area concerns the thermodynamics of 
the linear sigma- model. In ref. the finite-temperature version of the HFB-QRPA in 0(4) has been 
studied, in an attempt to describe the chiral phase transition. This transition was found to be of first- 
order. To leading order in the 1/A^-expansion, on the other hand, it turns out to be of second-order, 
as will discussed in future work. Similar to the 7T7r-scattering case it will be important to see what the 
next-to-leading-order corrections are. 
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7 Appendix 1 



In this appendix we want to show how one can handle safely the square roots originating from the 
Holstein-Primakoff mapping. Obviously one rather needs the expanded form of these functions in terms 
of the operators A and A + . Thus one could think of performing a formal Taylor expansion. However since 
these operators are non commuting one therefore must be careful when proceeding with the expansion. 
The major difficulty comes from the operator differentiation needed for the evaluation of each order in 
the Taylor expansion. Several methods have been designed for this purpose. Among these methods there 
exists a very simple and powerful one based on the parameter differentiation technique (see for instance 
|27|| ). Here we will present a slightly different derivation. 

The starting point is the parameter differentiation of the exponential. Given a function z of the complex 
variable A then the following identity holds: 

r l (96) 



dX Jo dX 



The identity above was derived by Snider in |2^] and can be recovered, following the same author, by 
first recalling the relation: 

~d\ ~ \ z ax z (97) 

m=0 

which can be shown to hold by a simple induction. Using this result the left hand side of eq.(p6[) reads: 

n=0 m=0 

Making now the substitution: p = n — m — 1, one observes that 

- — mlpl ^ = B(p+1, m+1) = ( 1 {l-a) m a p da (99) 
(p + m + 1)1 Jo 

where B(x,y) is nothing else but the beta function. Inserting now the above in eq. (|98|) one gets: 



<9e 2 ( A ) ^ ^ m\p\ z m dz z p 

d\ ^ n ^ (p + m+ 1)! ml dX pi 

m=0 p=0 KJr ' 

l z m dz 
ml d\ pi 



r 1 z m dz zP 

J ° n,p=0 



Now one can easily group the terms together and check that this is effectively the desired result of eq.(96). 
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Back to our problem of the square root expansion, we wish to make use of the identity derived so far. 
First to generate an exponential out of the square root one simply recalls the inverse Laplace transform: 



4V7TI Jc 



c+ioo 



e ts s 2 d s 



(101) 



Expanding the Holstein-Primakoff square root amounts then to expanding the exponential of the same 
argument times an arbitrary parameter s. For the purpose of the parameter differentiation we define the 
argument of the square root ( or of the exponential) as: 



t = t(X) = N + d 2 + \(dA + A + d + A + A) 



(102) 



For A = 1 one recovers the argument of the HP square root. We recall that the matrix d is diagonal, and 
that it does not commute with the operators A and A + . The fact that it scales like y/N is crucial 
because it helps to organize the expansion in the terms of powers of s/N which is the real parameter of 
the expansion (see text). 

Gathering everything together one can write finally the HP square root into the following Taylor expan- 
sion: 



t(A) 



j rc+ico 



q,p 



4^/7ri 



s 2 ds 



,at(A) 



A=0 



+ A 



Q e st(X) 



J A=0 



A 2 

+ Y 



Q2 e st(\) 



+ 



J A=0 



(103) 



As an application of the above considerations one can, through a trivial intermediate algebra, show that 
the lowest orders (used in this paper) for the HP square root read: 



(A+y/N + A+A) = N[d q ^l + d? q 5 q;p 



+ N2 



d 2 

u q ^ p 



X 1 ( d A + A+d) 

di V / q,p 



+ 



£ ^-r| d « { A+A ) q , p + Km (dA + AU) m4 



"m "'p 

d 



+ YF(q,m,p)^(dA + A + d) (dA + A + d) 

2 V / q,m \ / rn,p 



+ 



0(N~ 



(104) 



with the following definition of the function F{q,m,p): 



M+dl + Jl + dl-2 x /lTdt l Jl + dl-Jl + dl / ! ! 
^(q, m, p) = — — ,„ — ,„ — h ^ I — — + 



(di-dDidi-dj) 



d? q -d? p 



(d 2 q -cPJ (dj-cPJ^ 



(105) 
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One can further simplify these expressions by using the symmetry property of the bosons A, A + under 
the interchange of the indices q and p. The same expansion can be performed for the hermitian conjugate 
of the right hand side of eq.( |104] ). These two expansions are the building blocks of the Hamiltonian in 
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